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I— I, Abstract 

CN ■ A dark sector — a new non-Abelian gauge group Higgsed or confined near the GeV 

■ scale — can be spectacularly probed in low-energy e"'"e~ collisions. A low-mass dark 
i sector can explain the annual modulation signal reported by DAMA/LIBRA and the 

■ PAMELA, ATIC, and INTEGRAL observations by generating small mass splittings and 
• ■ new interactions for weak-scale dark matter. Some of these observations may be the first 

, signs of a low-mass dark sector that collider searches can definitively confirm. Produc- 

I tion and decay of 0(GeV)-mass dark states is mediated by a Higgsed Abelian gauge 

' boson that mixes kinetically with hypercharge. Existing data from BaBar, BELLE, 

^ . CLEO-c, and KLOE may contain thousands of striking dark-sector events with a high 

^ I multiplicity of leptons that reconstruct mass resonances and possibly displaced vertices. 

^ • We discuss the production and decay phenomenology of Higgsed and confined dark sec- 



tors and propose e'^e collider search strategies. We also use the DAMA/LIBRA signal 
to estimate the production cross-sections and decay lifetimes for dark-sector states. 
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1 Introduction 

Low-mass particles neutral under the Standard Model are nearly unconstrained by existing 
searches. We consider low-energy collider signatures of a new (9(0.1-10 GeV)-mass "dark" 
sector with gauge group Gd and light matter charged under it. The discovery of a dark sector 
that extends the structure of known low-energy gauge interactions would open a new frontier 
for particle physics and provide a second laboratory for fundamental questions including 
gauge unification and supersymmetry breaking. 

If Gd contains an Abelian subgroup U{l)n, high-scale physics generically induces kinetic 
mixing between the dark U{1)d and hypercharge, as sketched in Figure [l} this mixing medi- 
ates direct production of dark-sector matter in high-luminosity scattering experiments. The 
decay of the dark-sector matter can lead to spectacular events with high lepton and hadron 
multiplicities, lepton pair masses reconstructing resonances in the dark sector, and possi- 
bly displaced vertices. In optimistic scenarios, thousands of these events could be sitting 
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Figure 1: We consider a dark sector with non-Abelian gauge group Gu, which is Higgsed 
or confined at 0{MeV — 10 GeV). We assume that Gd contains a Higgsed Abehan factor 
U{1)d, so that the dark sector interacts with Standard Model matter through kinetic mixing 
of hypercharge with the U{l)r) gauge boson A' , of mass in the same range. Either the Higgsed 
or confined phases of Gd necessarily include new states that can be produced through A' 
interactions. 



undiscovered in data collected by BaBar, BELLE, CLEO-c and KLOE. Reconstructing new 
resonances in these events would reveal the dynamics of the dark sector. 

Evidence for a low- mass dark sector is emerging from a surprising source: accumulating 
hints from terrestrial and satellite dark matter experiments indicate that dark matter is not 
an afterthought of the Standard Model's hierarchy problem, but instead has rich dynamics of 
its own. The local electron/positron excesses reported by HEAT ^Lj, PAMELA [2j,|3j, PPB- 
BETS ATIC [5], and others |6l [7] are suggestive of weak-scale dark matter interacting 
with a new light boson, for which the U{1)d is a natural candidate if its mass is 0{GeV). 
Dark-sector interactions can also generate the mass splittings among dark-matter states 
suggested by other experiments [HI [9]. Dark matter scattering inelastically into an excited 
state split by 0(100 keV) can simultaneously explain the annual modulation signal reported 
by DAMA/Nal [10] and DAMA/LIBRA [II] and the null results of other direct-detection 
experiments [121 [13]. Likewise, the INTEGRAL 511 keV excess at the galactic center appears 
consistent with the excitation of dark matter states, but requires a slightly larger splitting of 
0{MeV) [HIE]. 

The electron/positron excesses and the splittings suggested by the DAMA and INTE- 
GRAL signals independently motivate an (9(GeV)-mass dark sector. If any of these anoma- 
lies are signals of dark matter, the new dynamics required to explain them can be discovered 
at e+e^ colliders. Among the anomalies, DAMA's signal offers the most precise predictions 
for e~^e~ collider physics: the scattering rate is sensitive to the strength of kinetic mixing 
between the Standard Model and the dark sector, and gives reason to expect an observable 
direct production cross-section for the dark sector. 

Outline 

In the remainder of this introduction, we further develop the motivation for a kinetically 
mixed light dark sector, and briefly describe the resulting events in low-energy e"'"e~ collisions. 
In Section 11.11 we discuss the kinetic mixing that couples the Standard Model to the dark 
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sector. In Section ll.2[ we summarize the evidence that dark matter, with a mass of 100 
GeV to 1 TeV, may interact with a low-mass dark sector. In particular, we introduce two 
frameworks for inelastic dark matter that motivate the Higgsed and confined dark sectors 
considered here. Though the dark matter in each model is heavy, in either case the generation 
of inelastic splittings requires light matter in the dark sector that can be produced at e^e^ 
colliders. We describe the general structure of dark-sector production events in Section 11.31 

We then work our way into the dark sector, and back out. Two production modes are of 
interest, and their cross-sections are nearly independent of the detailed structure of the dark 
sector. We present the production cross-sections for general parameters in Section 12.11 and 
for the parameter ranges motivated by the inelastic dark matter models in Section 12.21 In 
Section [31 we discuss the decay modes of metastable states, and the shapes of events expected 
within either Higgsed or confined dark sectors. 

Section H] contains a brief discussion of search regions and strategies, which synthesizes 
the experimental consequences of the results of Section [HJ Section H] is self-contained, and 
readers who are not interested in the detailed decay phenomenology of different dark-sector 
models can skip Section [3], and only read this section. We conclude in Section [51 

Appendix [X] reviews the existing constraints on a new U{1)d gauge boson with mass 
0{GeV) kinetically mixing with hypercharge. All numerical results are presented both for 
arbitrarily chosen fixed couplings, and for couplings normalized to the DAMA/LIBRA inelas- 
tic dark matter scattering cross-section using the prescription in Section 12.21 and Appendix 



The B-factory phenomenology of pure U{1) dark sectors has been discussed in [16], which 
notes a multi-lepton final state from a "Higgs'-strahlung" process, and the phenomenology 
of MeV-scale dark-matter models in which 7 + e~^e~ or 7 + nothing final states dominate 
has been considered in [T7]. The discussion of Section 3 is similar in spirit to the hadron 
colhder phenomenology for Hidden Valleys [THl [T^l [201 [2T] ^md the Higgsed dark-sector models 
discussed in [H [221 [23] (see also [21]). Other work that discusses new physics beyond the 
Standard Model leading to rare T decays at low-energy e~^e~ colliders (albeit with very low- 
multiplicity or even no leptons in the final state), such as T — *■ 7 -(- nothing, T —>■ invisible, 
or T — >^ 'jH (where if is a light Higgs boson) , includes [231 ESI [271 EB] ■ 

1.1 ^(1) Mixing with Dark Sectors 

Any new Abelian gauge group U (1) n has a gauge-invariant kinetic mixing term with Standard 
Model hypercharge [291 [30] . 



Several high-scale mechanisms, summarized below, generate ey ~ 10~^ — 10~^ (we note that 
from string theory the range for ey can be even larger [3T]). Therefore, the setup illustrated 
in Figure [H is generic if there is any new gauge group containing a t/(l)_D factor. 

The operator ([1]) is not generated at the Planck scale if f/(l)y is embedded in a Grand 
Unified Theory (GUT), because it is not SU{5) gauge-invariant. However, Planck-suppressed 
gauge-invariant operators involving GUT-breaking Higgses generate mixing 
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when these Higgses obtain vacuum expectation values, where the power p depends on the 
representation content of the Higgses. For example, p = 1 for the operator -^Tr[$F^^]F^'^, 
where $ is an SU{b) adjoint Higgs, and and are the 5'f/(5) and dark sector field 
strengths, respectively. 

Kinetic mixing is also generated by heavy split multiplets charged under both U{l)j:, 
and S'f/(5). Integrating out these multiplets generates loop-level mixing proportional to the 
logarithm of mass ratios within the multiplet: 

9i9d m \ ( ^\ 
Qoo.~^iV,log^-J. (3) 

The mass splittings can themselves be generated by either tree- or loop-level effects, leading 
to ey suppressed by one or two loop factors. 

At low energies, the kinetic mixing can be removed by a field redefinition of A^'^, inducing 
e-suppressed couplings of the A to electromagnetically charged states, 

C D egnA'^J^M, (4) 

where e = eycos^vy and Oyy is the Weinberg weak mixing angle. Therefore, a light f/(l)_D 
gauge boson can and naturally does open a path into and back out of any dark sector. 

1.2 Dark Matter Motivation for a Light Sector 

Searches for a GeV-mass dark sector are particularly motivated by recent evidence for 
electroweak-mass dark matter with non-trivial structure. The consistency of several inde- 
pendent observations of distinct phenomena with the same qualitative framework is striking. 
Moreover, if any one of these observations is found to have a non-dark-matter origin, a 
low-mass dark sector is still strongly motivated by the other data. 

Satellite observations and terrestrial experiments provide two kinds of evidence for a 
low-mass dark sector: 

1. Positron and/or electron flux measurements by PAMELA [21 E], PPB-BETS |;4j and 
ATIC [5] point to an unknown local source of high-energy (100 GeV-TeV) electrons and 
positrons. If their source is dark matter annihilation or decay, synchrotron radiation 
from these electrons and positrons could also explain the WMAP haze near the galactic 
center [32j. Two features of these experiments are suggestive of light states coupled to 
dark matter: 

• The excess of electrons and positrons, without a visible anti-proton excess |3], 
suggests dominantly leptonic decay or annihilation channels. One way to assure 
this is if the dark matter decays or annihilates into 0{GeV) particles that are 
kinematically forbidden from decaying to baryons [8]. 

• If the excesses are produced by dark matter annihilation, the annihilation cross- 
section required to explain the signal is 10-1000 times larger than the thermal 
freeze-out cross-section. This can be explained by Sommerfeld enhancement of 
the annihilation cross-section through an 0{GeV) mass force mediator at low 
velocities [S]. 
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2. As we will discuss below, non-Abelian gauge interactions Higgsed or confined at the 
GeV scale can generate small mass splittings among states charged under them. Two 
observations suggest that dark matter may have such splittings: 

• The INTEGRAL telescope [H] has reported a 511 keV photon signal near the 
galactic center, indicating a new source of ~ 1-10 MeV electrons and positrons. 
This excess could be explained by annihilation of light (9(1-10 MeV) dark matter 
[33], or by C(100 GeV-1 TeV) dark matter with O(MeV) excited states [15]. In 
the latter case, dark matter excited by scattering decays back to the ground state 
by emitting a soft e^e~ pair. 

• The DAMA/Nal [lO] and DAMA/LIBRA [TT] experiments have reported an an- 
nual modulation signal over nearly eleven years of operation. Modulation is ex- 
pected because the Earth's velocity with respect to the dark matter halo varies 
as the Earth moves around the sun, and the phase of the observed modulation is 
consistent with this origin. A simple hypothesis that explains the spectrum and 
magnitude of the signal, and reconciles it with the null results of other experiments, 
is that dark matter-nucleus scattering is dominated by an inelastic process, 

XN^X* N, (5) 

in which the dark matter x scatters off a nucleus into an excited state x* with 
mass splitting 5 ^ 100 keV [T2] . 

Among these hints of new dark matter interactions, the inelastic dark matter (iDM) 
interpretation [12] of the DAMA/LIBRA signal deserves particular attention, because it is 
sensitive to the strength of coupling between the dark sector and the Standard Model. While 
the other data discussed above is consistent with extremely small kinetic mixing e, the rate 
of scattering at DAMA/LIBRA fixes e within a factor of few for fixed model parameters 
and A' mass, in a typical range 10~^ ^ ^ ^ 10~^. Therefore, the spectrum and rate of the 
signal reported by DAMA/LIBRA can be translated into predictions for the couplings of 
light states observable at e^e~ colliders, as we will discuss below and in Section 12.21 

We summarize here how iDM reconciles the DAMA/LIBRA signal with null results from 
other searches. Inelastic scattering can only occur when the center-of-mass kinetic energy of 
the dark-matter-nucleus system exceeds the splitting 5. This implies that for a given target 
nucleus mass m^v and recoil energy Er, the minimum lab- frame velocity of dark matter 
required for inelastic scattering is 

fc„(£„)=.^=l=(!^ii^ + ^), (6) 

V2m7v^i? V / 

where /i is the dark-matter-nucleus reduced mass. This energy threshold enhances the sen- 
sitivity of the DAMA/LIBRA experiment over other direct detection experiments in three 
ways: 

1. Dark matter with velocity below v^am — a/2(5//x does not scatter. This minimum 
velocity increases for small m^r, thereby enhancing scattering off heavy nuclei such 
as Iodine used in DAMA/LIBRA relative to scattering off Silicon and Germanium at 
CDMS [M]. 
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2. Dark matter direct searches obtain the tightest bounds on elastically scattering dark 
matter from low-energy nuclear recoils. If dark matter scatters inelastically, low-energy 
recoils require higher dark matter velocity, and are suppressed — see equation ([6]). 

3. The fraction of the halo velocity profile that exceeds fmin — a/25/ fi is exponentially 
sensitive to changes in the Earth's velocity, leading to an 0{1) annual modulation of 
the scattering rate, in contrast to the ~ 2 — 3% modulation typical of elastic scattering, 
which arises from the simple linear dependence of the dark matter flux on the Earth's 
velocity. 

In this paper, we will consider dark-sector gauge groups Higgsed or confined at the GeV 
scale, both of which can generate (9(100 keV) splittings among (9(100 GeV-TeV) dark matter 
states, with a kinetically mixed U{l)r) mediating inelastic scattering off nuclei [SI E]- If the 
dark gauge group Go contains a Higgsed non-Abelian factor, radiative effects can split all 
components of matter charged under Gn, with splittings 

6 ~ aD^mw^ (Higgsed/radiative), (7) 

where Amwo is a splitting of dark-sector gauge boson masses [8]. For example, Amwo ~ 
~ 1 GeV and an ^ 10~^ gives 6 ~ 100 keV. 
If instead a non-Abelian factor of Gn confines at a scale A^, a heavy- flavor meson can 
be cosmologically long-lived and thus a dark matter candidate [5] . A particularly interesting 
limit for inelastic dark matter is when one dark quark is heavy (mg^ ^ ttldm ~ (100 GeV — 
1 TeV) ^ A^)), and the other light {rriq^ ^ -^d)- A cosmologically long-lived Qd forms 
uncolored hadrons after confinement, of which the lightest, a spin-0 QdQd meson, is a viable 
dark matter candidate. The spin-1 excited QdQd meson is split by hyperfine interactions, 
which give 

5 ~ — (confined/hyperfine), (8) 

where Nq is the number of colors of Gd. For example, ttidm ~ 1 TeV, Nq = 3, and A/) ~ 500 
MeV gives 5 ~ 100 keV. 

For both Higgsed and confined scenarios, the f/(l)D described above in Section [LT] can 
mediate inelastic scattering of dark matter off Standard Model nuclei, with elastic scattering 
naturally suppressed. These mechanisms can also generate a splitting of 0{MeV) to explain 
the INTEGRAL excess [IS]. 

Remarkably, though the mechanisms for generating the splitting from Higgsed and con- 
fined dark sectors are quite different, both are suggestive of dark sectors with multiple new 
states at (9 (GeV) — precisely the range accessible in a variety of low energy e~^e~ machines. 
In the Higgsed case, these states include gauge and Higgs bosons. In the confined case, these 
states include glueballs, light-flavor mesons, and baryons. Therefore, if either scenario is cor- 
rect, experiments such as BaBar, BELLE, CLEO-c, and KLOE may be capable of discovering 
the kinds of spectacular events that we consider here. 

We will see that the combination a^e^/m^,, where = 5'£)/47r, determines the DAMA/ 
LIBRA modulation rate. Therefore, we can use the observed modulation signal to estimate 
a£,e^ as a function of m^/, which in turn determines the expected production cross-section 
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Figure 2: Left: Cartoon of an event in which non-Abehan gauge bosons in a Higgsed dark 
sector are produced through an off-shell A: the gauge bosons may decay into other dark 
bosons (or fermions if present), which in turn decay to light Standard Model fermions. Right: 
a similar final state recoils off a hard photon in A' + 7 radiative return production. 

for ?7(l)D-charged matter at e^e~ colliders. We discuss this estimate in Section [2.21 The 
expected cross-sections are (9(fb) for m^/ ~ GeV. For heavier A', larger mixing parameters are 
required to reproduce the DAMA/LIBRA scattering rate; the resulting expected production 
cross-sections scale as ~ m\i. Therefore, very large cross-sections are expected for a few-GeV 
A', while a dark sector with sub-GeV N mass may evade detection in existing data. 

1.3 Striking Events at B-factories 

Existing B-factory datasets, containing over 1.4 ab"^ of data altogether, are ideally suited to 
search for a dark f/(l)D in the ~ 100 MeV - 10 GeV mass range, and for any dark sectors 
they connect to. There are two leading production processes: an A! with mass beneath 10 
GeV can be produced on-shell in association with a photon (the radiative return process). 
The A initiates a cascade of dark sector decays if kinematically allowed; otherwise it decays 
to Standard Model fermions or hadrons. In addition, any light charged states within the 
dark sector can be produced through an off-shell A carrying the full center-of-mass energy 
of the t^t~ collision. The signatures of these processes are quite spectacular because decays 
or parton showers within the dark sector generate high particle multiplicities. 

We focus on two cases, motivated by the dark matter models discussed in Section II. 2t 
a Higgsed dark sector containing only the non-Abelian gauge bosons and Higgses, and a 
confined sector with one light flavor. 

In a Higgsed dark sector, the fate of the initially produced dark states is determined 
by spectroscopy: they decay to lower-mass states within the dark sector if such decays are 



7 





SM final states 
in collimated jets 



^^ ^^^ \ > 



roughly spherical 
event 



boosted products 
recoil off photon 

v'd; 




Figure 3: Left: Cartoon of an event in which quarks in a simple confined dark sector are 
produced through an off-shell A'. The quarks shower and hadronize into mesons, which 
decay into Standard Model particles. Final states frequently contain many leptons, but can 
also include hadrons and long-lived dark states that escape the detector unobserved. Right: 
Phase space structure of different kinds of events. An off-shell A' produces jet-like structure 
if A/) <^ y/s (top), and approximately spherical final states if < y/s (middle). In A'/'j 
radiative return production, the dark-sector final state recoils against a hard photon. 

kinematically accessible and allowed by symmetry. These cascades eventually produce light 
states that can only decay to Standard Model final states. Gauge boson decays are suppressed 
by two powers of the mixing parameter e, and can be prompt or displaced. Higgs decays can 
be suppressed by e^, depending on kinematics, in which case they leave the detectors before 
decaying to visible matter. Typical events in a Higgsed dark sector can produce between 4 and 
12 Standard Model particles, with leptons being a significant fraction and easiest to observe. 
Caricatures of these events, with and without a recoiling photon, are shown in Figure El The 
decay phenomenology is similar in hadron colliders [21], [23] , and the pure Higgsed Abelian 
case for B-factories has been discussed in [16], though the dominant production modes and 
kinematics considered here are quite different. 

If the non-Abelian factor of Gd is confined, then the physical picture is very similar to 
the hidden valley models discussed in [THl [13 [SO]- U{1)d mixing mediates production of a 
light quark-antiquark pair in the dark sector. These states shower and hadronize, producing 
few dark-sector mesons with a roughly spherical distribution if the ratio y^/A/) is 0{1), 
and collimated jets if this ratio is large. Unlike the Higgsed scenario, the multiplicity of 
mesons in a typical final state is determined by the ratio of the production energy to Ad, 
not by spectroscopy. Different scenarios are caricatured in Figure [3l These sectors contain 
light mesons that can only decay to Standard Model final states. A single event can contain 
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a combination of prompt and displaced decays, and states that escape the detector. In 
particular, dark pions that cannot decay within the dark sector are very long-lived because 
Standard Model matter is neutral under U{1)d- This is quite different from the decays of 
hidden-valley pions, mediated by a Z' under which Standard Model matter and hidden- 
valley matter are both charged. Even when all states decay promptly, the reconstructed 
energy may still be less than the collider energy if some of the many tracks in these events 
are not reconstructed. Higgsed and confined-sector decays are discussed more carefully in 
Section [31 

These considerations motivate searching for dark-sector events in several channels at e+e^ 
colliders, as discussed in Section HI We will argue that six qualitatively different searches 
cover a wide range of possible phenomenology: 

1. M (exclusive), reconstructing Ecm (also 4£ + 7) 

2. M (exclusive), with displaced dilepton vertices (also 4£ + 7) 

3. > 5i + tracks (inclusive), reconstructing Ecm (also + 7) 

4. > 5i + tracks (inclusive), with displaced dilepton vertices (also + 7) 

5. Very high track multiplicity, with many tracks consistent with leptons 

6. 7 + nothing 

The combined results of searches in these regions can discover or exclude kinetically mixed 
dark sectors over a wide range of mass scales and couplings, including much of the parameter 
region motivated by inelastic dark matter. 



2 Production Modes and Cross-Sections at e^e Col- 
liders 

In this section, we discuss the production of dark-sector particles at e~^e~" colliders, with 
particular focus on the B-factories, BaBar and BELLE. The BaBar experiment at the SLAC 
National Accelerator Center collided e^e~ pairs, obtaining 430 fb~^ of integrated luminosity 
on the T(4S) resonance, at a center-of-mass energy of 10.58 GeV. The Belle experiment at the 
KEK laboratory in Japan obtained an integrated luminosity of 725 fb~^ on the T(4S). The 
two experiments also acquired ~ 270 fb~^ on the T(3S) and T(2S) resonances, and at nearby 
energies. The ~ 1.4 ab~^ of data available to the two B-factories make these experiments 
ideally suited for searching for the 0{GeV) dark sectors suggested by direct detection and 
astrophysical data. 

Over a large range of parameters, the cross-sections for the production of dark-sector 
particles scale as 

(9) 
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where i?cm is the center-of-mass energy of the colhder, ao = and go is the A' gauge 
couphng constant. The search sensitivity of a given e~^e~ machine above mass threshold 
scales as the ratio of integrated luminosity over squared center-of-mass energy, C\-at/El^. 
LEPI (Ecm ^ 91 GeV and ^int ~ 0.5 fb-^) and LEPII (Ecm ^ 189-209 GeV and Ant ~ 
2.6 fb~^) are much less sensitive to direct production of low-mass dark sectors than the B- 
factories. As we will review below, however, searches for rare decays of the Z-boson can 
probe some of the relevant parameter space. 

Compared to B-factories, lower-energy colliders such as DA$NE, which runs at E^.^ ~ 
m$ ~ 1.02 GeV, require much lower integrated luminosity to reach similar sensitivities to 
sub-GeV dark sectors. For example, the ~ 2.5 fb~^ of data collected by KLOE at DA$NE 
should make this experiment competitive with any B-factory searches for sub-GeV mass 
dark sectors, though ~ 1/3 as many dark-sector events would be expected at DA<I>NE than 
at BELLE. Importantly, the signatures of confined dark sector models can be quite different 
at E 

cm ~ 10 GeV compared to E^^ ~ 1 GeV. Since there is less showering of dark sector 
states for lower -Ecm? events at KLOE may be easier to reconstruct than B-factory events. 
We will nevertheless concentrate on the B-factory phenomenology in this work, with the 
understanding that at least part of the discussion is qualitatively unchanged for several other 
e+e~ colliders. 

In Section 12.11 we quantitatively discuss the production modes and cross-sections for 
general low-mass dark sectors coupled to the Standard Model sector via kinetic mixing. In 
Section [2721 we will assume that stable matter in the dark sector explains the direct detection 
signal reported by DAMA/LIBRA, within the framework of inelastic dark matter. Under 
this assumption, we will estimate the expected production rates at B-factories from the 
direct detection rate reported at DAMA/LIBRA. Theoretically reasonable parameter ranges 
permit cross-sections large enough to produce hundreds to tens of thousands of events at the 
B-factories. 

2.1 Production Processes for Low-Mass Dark Sectors 



At e^e~ colliders, there are two important dark-sector production modes. For definiteness, 
we focus here on the production of a pair of dark fermions, but the production cross-section 
for a pair of dark gauge bosons in a Higgsed dark sector, which occurs through mass mixing 
of the A' with non-Abelian gauge bosons, is similar. 

We first consider the direct production of a pair XX of dark-sector fields through an 
off-shell A', shown in Figure |l](a). The cross-section for this off-shell production process is 
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(10) 



Here, Nc is the number of colors in the dark sector gauge group Gd-, Nf is the number of 
dark sector particles Xi coupling to A' with charges qt, and the other variables were defined 
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Figure 4: Production modes of light dark-sector particles XX at B-factories. Left: Produc- 
tion through an off-shell A'. Right: Production of an on-shell A' and a photon — we assume 
the A' subsequently decays into lighter dark-sector particles. 



in equation (Q. This cross-section is non-zero even for rriA' much larger than Ecm, as long as 
mx < -Ecm/2. 

The second production mode, shown in Figure|4](b), is that of an on-shell A' and a photon, 
with differential cross-section 



2/1 4m^,/_B^^ 
9„^2„,2 / ™2 \ 1 + COS + 7" 2 / f2 ^2 



dcOsO V ^cm/ (l + COSe)(l-COSe) 

where 9 is the angle between the beam line and the photon momentum. This has singularities 
as the photon becomes coUinear with the initial-state electron or positron. The singularity 
is cut off by the electron mass, but if we demand that the photon be detected in the electro- 
magnetic calorimeter in the range cos6'min < cos^^ < cos6'max, the cross-section is 



where 



Q^Yo { + COs6'max)(l - COS gmin 
V (1 + COS 6'min) (1 - COS 6'niax 




(14) 



is an enhancement from the t-channel singularity (the numerical estimate is obtained from 
the BaBar calorimeter range in the center-of-mass frame, —0.92 < cos(^mm) ^ 0.87 [36j), and 
we assume m^, <C in equation f[T^ . The on-shell process (IT^ is enhanced by ~ ^0 
relative to the off-shell mode fllOl) . 

For mA' ~ 1 GeV, ao = a, and e = 10~^, the relevant cross-sections are C(fb) at the 
B-factories. With a total integrated luminosity of ~ 1.4 ab~^, this corresponds to several 
hundred to a few thousand dark-sector production events! In Figure O we show the inclusive 
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Production Cross-Section via Off -Shell A' Production Cross-Section for On-Shell A' + 7 
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Figure 5: Left: Inclusive cross-section at the B-factories (-Ecm = 10.58 GeV) for production 
of dark-sector states, X, via an off-shell A' as a function of ao^"^ and m^', for nix = 1 GeV, 
Nf = 1, Nc = 3, and unit charges qi (equation ffTOj) ). Note that the cross-section scales 
linearly with the number of dark flavors Nf and dark colors Nc- Right: Cross-section at 
the B-factories for production of an on-shell A' and a photon as a function of and m^/ 
(equation f[T^ ). Black lines correspond to fixed cross-sections. Also shown are the constraints 
on the couplings of a new U{1)d mixing with hypercharge from measurements of the muon 
anomalous magnetic dipole moment (shaded regions) [35]. The green dotted lines correspond 
to the lower bounds on the range of couplings that could be probed by a search at LEPI for 
rare Z-decays to various exotic final states, assuming that branching ratios as low as 10~^ 
can be probed. More details are given in Appendix \M 



cross-sections at the B-factories for off-shell and on-shell production of dark-sector particles 
as a function of either a£)e^ or and m^/. Figure [6] shows a similar plot for DA$NE. 

Figures E] and [H] also show constraints on the couplings of a new U{1) d mixing with hyper- 
charge. Shown in shaded regions are constraints derived from the muon anomalous magnetic 
dipole moment as in [35]. The measured value of is larger than the theoretically pre- 
dicted value, based on e"'"e~ annihilation to hadrons, by (302 ± 88) x 10^^^ [371 [38], and the 
constraints shown assume that the contribution from loops of the A' to is smaller than 
(302 + 5 X 88) X 10~^^. This constraint is discussed in more detail in Appendix \K\ 

Since any mixing between f/(l)_D and f/(l)y induces a coupling between Z-bosons and 
matter in the dark sector, decays Z XX are possible. Depending on the decay modes 
of X, which we will discuss in detail in Section [3], a variety of signatures are possible. For 
example, if the X's decay invisibly or outside the detector, the constraint is not better than 
the uncertainty on the width of the Z, which is about 0.1%. More exotic Z decays into 
multi-leptons are possible, but require dedicated searches. The best constraints on rare Z 
decay branching ratios are no better than the 10~^ — 10~^ level [SH], and even dedicated 
searches may not reach this level for complicated final states. In Figures [5] and [6l we show 
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Production Cross-Section via Off -Shell A' Production Cross-Section for On-Shell A' + 7 
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Figure 6: As in Figure El but for the center-of-mass energy of DA$NE (-Ecm = 1-02 GeV). 

(green dotted lines) the value of ao^'^ corresponding to a branching ratio of Z — s> XX as 
small as 10~^. The potential reach of rare Z-decay searches are discussed in more detail in 
Appendix \^ 

Finally, we mention two collider constraints on the production cross-section which are 
only relevant for certain decay modes. The first one is from the CLEO collaboration which 
searched for decays 'T(IS') ^ 7 + nothing' [IQ]. This is relevant for the case when the X's 
are long-lived and decay outside the detector. As we will see in Sections [3] and HI this is 
a very important search mode. Using an integrated luminosity of 48 pb~^ of data, CLEO 
found no peak in the photon energy in the range 1-4.7 GeV ^0]. The resulting limit can 
be interpreted as excluding cross-sections as low as 200-500 fb, depending on the A mass, 
though [IQ] assumes a different photon angular distribution than would result from dark- 
sector production. 

The second collider constraint comes from a search for T(3S') — > 7 A' — >■ 7/i"'"/i~ by the 
BaBar collaboration |4]j|, which is only relevant if A' decays to dark-sector particles are 
kinematically forbidden, in which case it decays to Standard Model states. In 30 fb~^ of 
data containing ~ 122 x 10^ T(3S') events, a 90% C.L. upper limit of roughly 4 x 10~^ on 
the 7/i~'"/i~ branching fraction was found for uia' of a few hundred MeV. This search would 
also be sensitive to about 500 7/i"'"/i~ events through v4'/7 production. In this mass range, 
Br{A' —y ~ 0.3 — 0.5 These results exclude production of about 1000 7^4' events, 

corresponding to a cross-section of about 30 fb. For A' masses of a few GeV, but below the 
r threshold, the limits exclude 250 7/2 events, but Br{A' —>■ fifi) ^ 0.2 — 0.25 so that a 
similar cross-section limit is obtained. We note that even if the A' is the lightest state of the 
dark sector, any other kinematically accessible states are pair-produced with a comparable 
rate through an off-shell A', which produces a signal in the multi-lepton channel (see Section 
E]), which has much lower QED backgrounds. Lower cross-sections can thus be probed by 
considering different final states. 
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In summary, the potential search mode for rare Z decays and the constraint do not 
exclude parameter regions in which thousands of events could be seen. Collider searches for 
'7 + nothing' and 7 + [i^ \x~ final states have probed some of this parameter space, but are 
quite model-dependent. More general final states could easily have been missed by existing 
searches, even if they occur with large cross-sections. 

Before discussing the decay of X's in Section [3l we use the hypothesis that an (9(100 
GeV) stable particle charged under Gb explains the direct detection signal reported by 
DAMA/LIBRA. We use this to estimate the dark sector couplings a^e^ and resulting pro- 
duction cross-sections of low-mass dark sector states at B-factories. 

2.2 Inelastic Dark Matter and Inclusive Cross-Sections 

Inelastic dark matter (iDM) can reconcile the annual modulation signal reported by the 
DAMA/LIBRA collaboration with the null results of other experiments [121113] . In proposals 
to generate the iDM splitting from new non-Abelian gauge interactions in a ~ GeV-mass dark 
sector, it is also natural to assume that scattering is mediated by a f/(l)D gauge boson N 
that mixes kinetically with the photon, with mass thai of (9(100 MeV - 10 GeV) [8]. In such 
models, the signal rate and spectrum reported by DAMA/LIBRA constrain the couplings 
of the A! to the Standard Model electromagnetic current. This allows us to estimate the 
inclusive cross-section for production of dark states at B-factories for a given iDM model. 

At energies beneath m^/, the dark matter interacts with charged matter through an 
effective coupling 

r ~ ^^^^ t n ^\ 

where J^^^Dark is a current composed of dark sector fields, and e and gB are the photon and A 
coupling strengths, respectively. The combination (az? = 9dI^t^^ controls the scattering 

rate of dark matter in direct detection experiments (see Appendix [B]), while determines 
the inclusive production cross-section at B-factories. Consequently, the production cross- 
section can be estimated as a function of m^/. 

We will consider two classes of iDM theories — one where the dark matter is directly 
charged under the IJ{\)d [H], and one where the dark matter is a neutral composite meson, 
built out of charged spin-1/2 constituents, one of which is light (~GeV) and the other heavy 
{jriBM ~ 100 — 1000 GeV) [9]. Higgsed dark sectors are an example of the first class of 
theories, in which the breaking of the dark gauge group generates mass splittings for the 
dark matter gauge eigenstates. The resulting mass eigenstates couple off-diagonally to the 
A', leading to inelastic scattering. The low-energy scattering matrix elements are 

(X* (pO I -^^Dark I X(P)) Charged Scattering = {V ^ VY , (16) 

where we have neglected higher order spin-dependent terms. Confined dark sectors are an 
example of the second class of theories, in which hyperfine interactions among the constituents 
split the ground state meson into a pseudo-scalar and vector If the J^^rk coupling breaks 
parity, then the leading scattering operators mediate inelastic hyperfine transitions (dipole 
scattering) with matrix elements 

(XV (P',e^) I ^Dark I X5(P)) Dipole Scattering ~ "T^ (P + P')'' ^ " 9> (17) 



14 



DAMA-Normalized ao e (Form-Factor iDM) 



DAMA-Normalized ao e (Charged iDM) 




1.0 2.0 

M^<GeV) 



1.0 2.0 

MA<GeV) 



Figure 7: odc^ as a function of m^/, normalized to the DAMA/LIBRA annual modulation 
signal. Black, red, and orange bands from top to bottom correspond to ttidm = 2000,800, 
and 200 GeV, respectively. Uncertainties in fitting to DAMA/LIBRA, as a function of the 
inelastic dark matter mass splitting, are reflected in the thickness of the bands. The left 
plot shows the normalized couplings for inelastic dark matter that scatters via form-factor- 
suppressed dipole scattering (equation (IT7|) ). typical for a confined dark sector, with C = 0.01 
(a£)e^ scales as 1/C*). The right plot shows the normalized couplings for charged scattering 
(equation f|T6l) ). which is typical for a Higgsed dark sector. 



where is the polarization of the outgoing xv and q = p—p' is the momentum transfer. The 
natural magnitude of the dipole moment b = ^ GeV~^ is the size A^^ of the composite 
dark matter [Ad is of order the confining scale of the dark gauge group). Relative to the 
charged scattering matrix element f|T6l) . this scattering is suppressed by a g-dependent form 
factor. 

For either current, the cross-section to scatter off of a nuclear target of mass rriN and 
charge Z with recoil energy Eji is. 



da SvrZ^aaDe^mAT 



1 (charged (vector-current)) 



^ '^■2 (form-factor (dipole)) 

2 

where a = 1^, f is the dark matter- nucleus relative velocity, and F{Ept) is a nuclear form 
factor. As discussed in Appendix [Bl we compute the scattering rate from equations f|T6|) and 
(ITTj) . fit the reported DAMA/LIBRA rate, and thereby obtain an estimate for the couplings 
a^e^ as a function of m^^. 

From equation (ITS!) , we see that the dipole scattering cross-section is a function of the 
momentum-dependent form factor ^^^^ ^ tt- To estimate the coupling a^e^ for a given 
niA', we must thus choose a value for Ad- Ad can in principle be determined by fitting 5 
to the DAMA/LIBRA data, since 6 is proportional to A'^j/itldm- However, the constant of 
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Figure 8: Inclusive cross-section at the B-factories for production of dark-sector states, X, 
through an off-shell A' (left) and for production of an on-shell A' with a photon (right) as a 
function of rriA', after normalizing to the observed DAMA/LIBRA modulation rate assuming 
form- factor-suppressed dipole inelastic dark matter scattering. We set rriDM = 1 TeV. In 
the left plot, we set nix = 1 GeV, Nf = 1, Nc = 3, qi = 1, and show various choices of 
C = {0.01, 0.1, 1.0}, while in the right plot, we show = {1, a, 10"^} with C = 0.01. Note 
that the off-shell production cross-section scales linearly with the number of dark flavors Nf 
and dark colors N^. Also shown on both plots are constraints from as well as the rare Z 
decay sensitivity region — see Appendix |A] for more details. 



proportionality is model-dependent. Following we introduce a scale Ad, such that 



S ^ (19) 
rriDM 



In general. 



T (strong coupling) , ^ 

I a'A^i (coulomb regime) 

where the first case corresponds to a heavy meson of a SU{Nc) confining at Ad, and the 
second case corresponds to a meson in a weak-coupling coulomb regime with coupling a'. 
Ad/ Ad is naturally a small number in either case. We absorb the uncertainty in both A^i 
and Cin in a single proportionahty factor 

C^l^l, (21) 

which is naturally small {C ~ 0.01 — 1). The dark matter scattering rate is now exactly fixed 
by the inelastic splitting 6 and C, allowing the determination of a^ie^ from the DAMA/LIBRA 
scattering rate as a function of rriA'- For the charged scattering ( fT6ll . there is no analogous 
unknown parameter, and we can solve for Q;£)e^ directly from ([Tl 
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Production Cross-Section for On- Shell A'+y (Charged iDM) 
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Figure 9: As in Figure [8], but here the production cross-sections are normahzed to the 
observed DAMA/LIBRA modulation rate assuming charged inelastic dark matter scattering. 



Figure [7| shows the DAMA/LIBRA-normalized couplings as a function of m^/ for charged 
and dipole scattering, and for different dark matter masses. In Figure [HI we show the B- 
factory cross-sections using dark sector couplings normalized for dipole scattering. In Figure 
[9l we show cross-sections normalized using charged scattering. For m^' > 100 MeV, B factory 
cross-sections are expected to easily exceed ~ 1 fb for composite iDM theories. For theories 
where the iDM scattering is charged-current dominated, ~ 1 fb cross-section can occur for 
m^/ > 1 GeV. In general, the estimated cross-sections scale as m\,. 

We note that the e+ and/or e" data from, for example, PAMELA, ATIC, and PPB-BETS 
cannot be used in the same way as the DAMA data to predict the cross-sections at low-energy 
e+e~ colliders. The reason is that the e"*" and/or e~ data does not directly constrain the A' 
coupling with the electromagnetic current. Since the dark matter annihilates into on-shell 
A's that decay with a lifetime that is largely unconstrained by the data, the mixing parameter 
e is also unconstrained. However, indirect constraints on the dark matter annihilation cross- 
section from, for example, 7-ray observations [111 HS] and CMB measurements [44j constrain 
the maximum allowable Sommerfeld enhancement, and currently disfavor very low A' masses 
(~ 100 MeV or less). 

3 Decays of Dark-Sector States into Standard Model 
Particles 

In this lengthy section, we discuss the decays of dark-sector states produced at e~^e~ colliders. 
Readers interested in the experimental consequences of these decays can skip to Section HI 
which is self-contained. 

The kinetic mixing that allows dark-sector states to be produced at e+e^ colliders also 
allows them to decay back into Standard Model states, mainly leptons and light hadrons. To 
demonstrate the range of phenomenology possible in both Higgsed and confined dark sectors. 
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we have considered several kinematic limits that exhibit different properties. We will also 
briefly consider the impact of Higgs mixing on decays. We organize our discussion around 
the following questions: 

1. What are the metastable states of the dark sector that can only decay through processes 
that involve kinetic mixing? 

2. What are the decay widths of these metastable states? 

3. What are the multiplicities and kinematic properties of these metastable states? 

We consider Higgsed dark sectors in Section 13.11 and confined dark sectors in Section 13. 2[ 
In both cases, decays of spin-1 states through A' mixing, suppressed by e^, are generically 
prompt, while other species can have displaced 3-body decays or invisible decays suppressed 
by e^. Therefore, different dark-sector spectra can give rise to very different phenomenology, 
with all dark-sector states decaying promptly into the standard model, all states so long-lived 
that they escape the detector, or a combination of prompt, displaced, and/or invisible decays. 
Light fermions in the dark sector lead to events with more invisible decay products, as we 
discuss in Section [331 

Through most of this discussion we consider only the decays generated by the A'/j 
kinetic mixing, but other interactions between the Standard Model and dark sector, involving 
either Standard Model Higgs bosons or new heavy fields coupled to both sectors, are also 
possible. Though these have no observable effect on direct production, they can dominate 
over the kinetic-mixing-induced decay modes that scale as e^. These interactions could be 
discovered through the observation of displaced decays that are not induced by kinetic mixing 
— specifically, a two-fermion decay of spin-0 states (Section I3.1.3P and any decay of the 
lightest dark-sector fermion (Section l3.3p . 

3.1 Events from a Higgsed Dark Sector 

We assume that the dark sector gauge group Gd has been completely Higgsed, giving mass 
to all dark-sector gauge bosons. We also assume that all custodial symmetries have been 
broken, so that couplings Wi Wj Wk, Wi Wj hk and Wi hi hj are generated between all Higgs 
and gauge bosons. The custodial symmetry breaking ensures that the A' mediates inelastic 
transitions between the dark matter and its excited states, which is needed to explain the 
DAMA/LIBRA and INTEGRAL signals [8l[23]. As all mass eigenstates are mixed, we will 
not distinguish between the Abelian gauge boson A' and non-Abelian gauge bosons, but refer 
to all as Wd- 

The high-energy collider phenomenology of these models has been discussed in [t23j, and 
many of the same features control low-energy collider signatures as well. The special case 
of a pure U{1) dark sector was recently discussed in [16j. In Section [3.1.H we determine 
the metastable states and their decays mediated through kinetic mixing. We discuss the 
final-state multiplicities in Section [3.1.21 In Section [3. 1.31 we comment on other interactions 
that induce faster decays than those mediated by kinetic mixing. 
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3.1.1 Metastable States and their Decays through Kinetic Mixing 

We begin with a somewhat careful discussion of the metastable states — those that can only 
decay through kinetic mixing to Standard Model states, or equivalently those whose lifetimes 
vanish in the e ^ limit. States that can decay to final states within the dark sector are not 
considered metastable, and these decays are typically prompt. The Standard Model states 
produced in the decay of metastable states include mainly leptons and pions. We denote the 
mass of the lightest dark-sector gauge boson by min^ mass of the lightest Higgs 

by ^/lomin- Higgs heavier than 2mvy^^.^ can decay to two gauge bosons, and any Wd 
heavier than 2m/i^ ^.^^ can decay to two Higgses. Therefore, if mtj^ ^.^ > ^niwo only Wd^s 
are metastable; if min > ^m^^^.^, then only Higgses are metastable. More than one 
dark gauge or Higgs boson may be metastable in each case, but they will decay similarly. If 
^mwa min < min < '^''^Wd min 5 thcu thc dark sector contains at least one metastable Higgs, 
and at least one metastable gauge boson. 

— Metastable States and Their Decays 

Wn prompt decay to // (2-body 

displaced decay to ffWn (3-body e'^) 
[Wd prompt decay to // (2-body e^) 

'hu invisible decay to 4/ (4-body e^) 
Wd prompt decay to // (2-body e^) 

hu invisible decay to 4/ (4-body e*) 



2- 

1- 

1 __ 

2 



Figure 10: Summary of decay phenomenology in Higgsed dark sectors. While decays within 
the dark sector are prompt, the lightest dark-sector state can only decay into Standard Model 
leptons or hadrons with a width that is suppressed by some power of e. Generically, this 
decay is prompt if the lightest state is a gauge boson, while it is displaced or invisible if it is 
a Higgs. The figure summarizes the decay phenomenology as a function of the mass ratio of 
the lightest Higgs boson to the lightest gauge boson. 



The importance of this characterization is that the gauge and Higgs bosons have very 
different lifetimes: gauge bosons typically decay promptly, while Higgs decays give rise to 
displaced decays or escape the detector before decaying. A typical dark-sector production 
event will eventually decay to a collection of metastable states, consisting either of only Higgs 
bosons, only gauge bosons, or a combination of both. This leads to a variety of scenarios with 
very different decay phenomenology, summarized in Figure [TUJ In the following, we justify 
the decay properties claimed in that figure. 

If Wd^s are metastable, they can decay to Standard Model leptons and pions through 
kinetic mixing. By assumption, all gauge bosons have some mixing angle 6 with U{1)d, 
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giving rise to a decay width 




2m} ' 


)iVeff 




f 


Anij 





2.4eV(^) ^Mt^) Wl-^(l + ^)iV,ff, (22) 

where mi denotes the mass of a lepton or a hadron. Here A^cfT counts the number of available 
decay products: N^^ = 1 for mw^ ^ 2m^ when only Wd — > e"'"e~ decays are possible, 
and 2 + i?(mi4/Q) for m^o > 2m^, where i? is defined to be the energy dependent ratio 
a{e~^e~ — >• hadrons)/(j(e~*"e^ /W^/U ) [39j. 

The decay length of Wd in its rest frame is given by 

CT{Wn ^ ~ 8 X 10-6 cm ( ^] ' r ^ f ( , (23) 



and is thus prompt for 6 ~ The combination a^e^ can also be normalized to its 

DAMA/LIBRA expected value. In this case, 

ct{Wd ^ tr) ~ 10-^ cm 0-^ — h 24 



where m^M is the mass of the dark matter (see Section [2^ . We see that prompt decays are 
also generic. However, due the strong dependence of cr on l/m\,, displaced vertices are also 
possible for m^' ~ 0(100 MeV). 

If Higgses are metastable, they can decay to Standard Model leptons and hadrons through 
one or two off-shell Wd^s (heavier Higgses can have 3-body decays to a lighter Higgs and 
Wl) If m/i^ > mwominy final-state Wd can be on-shell, and the lifetime is 

2 

TihD ^ W,W; ^ W,tn ~ ^^mhjsimwjm^^), (25) 

where we have assumed m\Y^ = my/_. for simplicity. This has the same form as the h ZZ* 
partial width in the Standard Model. The function /3(a;), | < x < 1, is a 3-body phase 
space integral that can be found e.g. on page 30 in [IS]. The decay length in the rest frame 
of ho is shown in Figure [TT] for ao = ol and for both e = 10^^ and 10^*^. From the figure, 
we see that hn decays can be prompt, displaced or invisible, depending on the parameters. 
Similarly diverse decay scenarios are also shown in Figure [T2| where we plot the decay 
length after normalizing to its DAMA/LIBRA expected value. 

If rnhj^^^^ < "rnvKomin; then the lightest Higgs can only decay through two off-shell gauge 
bosons, which is dominated by a four-fermion final state, and is parametrically given by 

V{hn - W*r>W*r> - U) ~ (26) 

izvr 

A 2-body decay to two Standard Model fermions through a one-loop diagram mediated by 
I^d's is parametrically of the same size. Because of the strong suppression, such a light 
Higgs is expected to escape the detector, even for large an- Note that in this case all but the 
very lightest Higgs have three-body decays to either a lighter gauge boson or a lighter Higgs, 
so that it is not hard to construct dark sectors in which prompt, displaced, and invisible 
decays all occur. 
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Figure 11: Decay length (cr) of a dark-sector Higgs boson, ho — ^ WiW* Wii£, as a 
function of its mass m/^^ and the dark-sector gauge boson mass mwa: for an = « and 
e = 10~^ (left) and e = 10~^ (right). Here i denotes a Standard Model lepton or hadron. 
The light to dark (light-blue/blue/dark-blue) regions indicate decays that are prompt (< 50 
/xm), short displaced (< 3 cm), and long displaced (> 3 cm), while in the white region the 
decays are invisible (> 50 cm). 



3.1.2 Event Shapes in a Higgsed Dark Sector 

Specific choices for the structure of a Higgsed dark sector can dramatically change the typical 
decay chains within the sector. It is typical for dark-sector production to result in a high- 
multiplicity and lepton-rich final state. Here, we discuss two examples. 

A four-lepton final state is expected from production of a pair of metastable Wd bosons 
that decay directly to Standard Model fermions (Figure [T3r a)). The lepton kinematics de- 
pends on the Wn masses: if ^ V^; then the observed leptons will lie in two approxi- 
mately coUinear pairs. 

It is also possible to produce heavier gauge bosons that decay into a lighter W'j^ and ho of 
the dark sector, as in Figure [TBT b). One limit in which this is guaranteed is a dark sector with 
the structure of a charge-breaking Standard Model, with SU{2)d x U^IYj^ Higgsed to U{1)d 
at aiid with charge {U{1)d) being broken at a lower mass scale rriA'- Such models have 
an approximate global custodial symmetry f/(l) that suppresses the charged Wd's mixing 
with the A', so that the heavy bosons are produced far more efficiently than the lighter 
A' ^23j. In this case, a large number of leptons (and possibly pions) are expected. Likewise, 
it is possible to produce dark Higgs bosons directly (a special case is the "Higgs'-strahlung" 
discussed in [T6]). 
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Figure 12: Decay length (cr) of a dark-sector Higgs boson, ho WiW* Wii£, as a 
function of its mass m/^^ and the dark-sector gauge boson mass mvv^,, for an = <y- Here the 
combination a^e^ has been normahzed to its DAMA/LIBRA expected value, assuming that 
the dark matter (with mass 1 TeV) is part of the Higgsed dark sector and charged under the 
U{1)d (see Section [2.21 for further explanation). The shaded regions in this figure are the 
same as in Figure [TTl 



3.1.3 Decay Modes from Additional Mixing Interactions 

Additional couplings between the Standard Model and dark sector are possible, but they in- 
volve massive Standard Model-charged states and so are irrelevant to dark-sector production. 
They can, however, compete with the e^-suppressed decays of any light Higgs bosons in the 
dark sector. For example, the mixed Higgs quartic 

5C = K\hD\'\hsM\' (27) 

leads to Higgs decay through mixing with the Standard Model Higgs. If the theory is su- 
persymmetric at high energies, ~ ^ggo is generated by D-term mixing More 
generally, loops of link fields charged under Gd and SU{2)l x f/(l)y and loops of ordinary 
Standard Model matter will generate a quartic Higgs mixing through their Yukawa interac- 
tions with both Higgses. A rough upper bound on the mixing, 

A. < (28) 

can be obtained by demanding that rrth^ receives at most 0{1) corrections from this quartic 
term. If the U{1)d symmetry-breaking is generated by ( 1271) . then this limit is saturated. 
When ho gets a vacuum expectation value, the quartic coupling ( 1271) generates mixing 
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Figure 13: Two examples of cascade decay chains for Higgsed dark sectors. Left: The 
production of a pair of dark-sector gauge bosons Wd can lead to a four-lepton event. Right: 
The production of a pair of heavy dark-sector gauge bosons, which decay within the dark- 
sector to lighter gauge and Higgs bosons, can lead to an event with a large number of leptons 
(and possibly pions). 



between ho and hsM with angle 

e ~ - -^./f , (29) 

■^SM VsM ^SM V 

where the latter limit corresponds to saturation of the bound (128|1 . This mixing leads to a 
decay width 



nhr> - £+r)_ ~ y^mnj' - f m,,^ ( ^ ) , (30) 



where the final expression is again obtained by saturating ( !28|) . Assuming saturation, the 
decay length oi ho in its rest frame is given by 

where Ut ~ 0.01 is the r Yukawa coupling. This decay can thus lead to observable decays with 
C(cm) displacements if hr, is above the r threshold, and saturates the naturalness bound. 
It dominates over the three-body decay (l25l) if e < 10~^. Below the r threshold, however, it 
is unlikely to produce observable decays. 
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3.2 Events from a Confined Dark Sector 



In this section, we discuss the events obtained with confined dark sectors. We begin in Section 
13.2.11 with a review of the spectroscopy in confined dark sectors containing one or multiple 
light flavors, and discuss the shape of production events. As for the Higgsed case, the events 
can usefully be characterized by discussing the metastable states, i.e. those states that do 
not decay to other dark-sector particles. For multiple light flavors, these metastable states 
are either the dark-sector pions (Section [3.2.21) or the A' (Section [3.2.31) . depending on which 
is lighter. For a single flavor, the metastable states are either the A' (also Section [3.2.31) . or 
the tower of mesons lighter than 2m^^ , where m^^ is the mass of the lightest meson (Section 
I3.2.4p . Likewise, the lightest baryon is exactly stable and a tower of baryons lighter than 
+ ''^v'd metastable. 

3.2.1 Preliminaries: Spectroscopy and Event Shapes 

We begin by reviewing the spectra of mesons and baryons in a one-flavor model. A model 
with only one light flavor has no light pion; instead, the lightest meson is expected to be 
a pseudoscalar meson {J^^ = 0~~^), which we will call r/^ in analogy with the Standard 
Model 7]'. The second-lightest meson, again named by analogy with the Standard Model, is 
the ujd (1 )• Excited states and higher-spin mesons are also expected, but their spectrum 
is not known precisely. We will only ever be interested in the tower of mesons with mass 
less than '2m^i^. Mesons above this threshold can decay strongly to lighter mesons, and are 
quite broad. The one-flavor dark sector also contains a tower of baryons, beginning from the 
spin-A^c/2 A/), where Nc is the number of colors of the conflning gauge group. For small N^, 
these have comparable mass to the tower of light mesons. We expect that any metastable 
glueball states will mix with the mesons, and decay through this mixing. 

In a model with Nf > 1 light flavors, conflnement breaks an approximate SU{Nf)L x 
SU{Nf)ji X f/(l) flavor symmetry to a diagonal subgroup SU{Nf)ii x f/(l). There are then 
Nj — 1 light pions. All other mesons and (for Nc > 2) all baryons are parametrically heavier 
than the pions, and are produced much more rarely; therefore, we will focus on the pions. 

The flnal-state kinematics and multiplicities can be estimated from those in the Standard 
Model by analogy with 'e"'"e~ —>■ hadrons' processes and a simple rescaling of the QCD to 
dark-sector conflning scale, or inferred from parametrizations such as the quark production 
rule [38]. These properties depend dominantly on the ratio of the conflning scale to the 
energy y/sh going into the conflned sector {^/Sh = -Ecm for e~^e~ hadrons or e'^e~ 
dark-sector hadrons, while ^/sh = for the radiative return process). We will parametrize 
this by Re = \fshlmr(^. Thus, for example, from low-energy SPEAR data we see that jet 
structure is not observable for Re ~ 4 (3.6 GeV collisions), but is observable by Re ~ 8 jH]. 

Likewise, dark-hadron multiplicities can be estimated from QCD track multiplicities in 
e^e~ collisions; what is measured is the charged-track multiplicities inch)-, not the total 
hadron multiplicity, but studies of the total energy of reconstructed tracks suggest that ap- 
proximately half of produced hadrons are charged |49j . Doubling observed track multiplicities 
from e.g. [501 [51], 'we see that over the range of interest, hadron multiplicities range between 
(A^) ^ 10 at 7.6 GeV {Re ^ 8) to ^ 40 at 90 GeV {Re ~ 100). 

At low Re in a one-flavor model in which 77^ is the lightest state, we expect (A^) ~ 
^/s/nirii , so that the production of non-relativistic hadrons dominates, while the produc- 
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Figure 14: Decay length (cr) of the dark-sector scalar meson 77^, decaying through the axial 
anomaly to two on- or off-shell A' gauge bosons, which in turn decay into two Standard 
Model leptons or hadrons. Here an = d, and e = 10~^ (left) or e = 10~^ (right). The shaded 
regions in this figure are the same as in Figure [TTl This figure is also applicable to decays of 
light pions for dark sectors with multiple light flavors. 

tion of a much lower multiplicity of hadrons is exponentially suppressed (see e.g. [22] for 
experimental evidence of this effect in QCD events with very few pions). It is difficult to 
make more quantitative statements in this region of parameter space based on either direct 
extrapolation from QCD or parametrized models. The QCD extrapolation would suggest 
hadron multiplicities larger than Ecm/^ri'^-, which are consistent in theories with light pions 
but unphysical if the rj'j^ is the lightest state. Parametrizations such as the quark produc- 
tion rule [IB] can in principle be modified for the one-flavor scenario, but they deviate from 
experimental data precisely at small Re- 

3.2.2 Dark Sectors with Metastable Pions 

In a model with two or more flavors and a heavy A\ the lightest dark-sector states are pions. 
Since they are then kinematically forbidden to decay into any other dark-sector states, they 
are metastable. 

Kinetic mixing does not generate the most general decays allowed for the dark pion in an 
effective field theory. In particular, the helicity-suppressed decay of a "f-pion", vr^ — > 
which dominates in Hidden Valley models with a Z' that couples directly to both dark-sector 
and Standard Model matter, is not generated by kinetic mixing. This can be most readily 
seen before performing the field redefinitions that diagonalize the gauge field kinetic terms. 
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Figure 15: Decay length (cr) of the dark-sector scalar meson 77^, decaying through the axial 
anomaly to two on- or off-shell A' gauge bosons, which in turn decay into two Standard 
Model leptons or hadrons. Here = a, and the combination ane^ has been normalized to 
its DAMA/LIBRA expected value, assuming that the dark matter is a 1 TeV dark meson 
neutral under U{1)d, but consists of a light and a heavy quark charged under U{1)d (see 
Section [22]) • The shaded regions in this figure are the same as in Figure [TTl This figure is 
also applicable to decays of light pions for dark sectors with multiple light flavors. 

Instead, let us treat the kinetic mixing as an insertion 

A'jj Aem 

By Lorentz invariance, the pion current must have the form 

(0|A;|7r(p))cx^p^, (33) 

which vanishes when dotted into the kinetic mixing tensor. In other words, kinetic mixing 
only mixes the transverse modes of the two gauge bosons, to which the pion does not couple. 
The vanishing of this decay mode is obscured when the couplings are diagonalized. The A' 
couples to the photon current, Jem, but this vector-like current vanishes when contracted 
with p^. The A' also couples to the chiral current, Jz, but this contribution is precisely 
canceled by the Z's coupling to the dark photon current, Ja', so this also does not mediate 
decays. 

Instead, the pion can only decay through the anomaly diagram, but with one or both A' 
off-shell. The resulting decay width scales as 



26 




Figure 16: Decay length (cr) of the dark-sector vector meson uj£,^ which mixes with the A' 
and decays to two Standard Model leptons or hadrons. Here ao = en, and e = 10~^ (left) or 
e = 10~^ (right). The shaded regions in this figure are the same as in Figure [TTl 

where p = 1 for three-body decays and p = 2 for four-body decays, and the pion decay 
constant /^^ can be estimated by scaling the pion decay constant in the standard model, 
fiTo = fT^^r/j^/^-n'sM- These decay lifetimes are identical to those of the in a one- flavor 
model, and the rest frame displacements are shown in Figure [TH If uia' < itLt^j^ < 2mA', a 
three-body decay can give rise to displaced vertices, but this mass condition is not generic. 



When TTiA' > TT^-Ko^ the four-body mode is suppressed by e'^, leading to very large displace- 
ments: 



A loop decay to two leptons is also possible, but comparable in rate to the four-body decay (in 
fact, it receives chiral suppression ~ (mf/m,ro)^ relative to the four-body decay). Therefore, 
in this region of parameter space the pion escapes the detector without decaying, and is 
invisible. 

Assuming that dark-sector interactions break the dark flavor symmetries, a heavier pion 
that cannot decay to two lighter pions may still decay via an off-shell A' into a lighter pion 
and Standard Model leptons/hadrons. This decay is suppressed by a^e^ , and so is typically 
prompt or displaced. 

3.2.3 Dark Sectors with a Metastable Light A' 

In this section, we discuss the case when the A' is the lightest state in the dark sector. We 
consider both a one-flavor and a multi-flavor model. In both cases, final states consisting of 
prompt Standard Model leptons and hadrons are likely. 
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Figure 17: Decay length (cr) of the dark-sector vector meson uoo-, which mixes with the A' 
and decays to two Standard Model leptons or hadrons. Here ao = a, and the combination 
a^e^ has been normalized to its DAMA/LIBRA expected value, assuming that the dark 
matter is a 1 TeV dark meson neutral under U{1)di but consists of a light and a heavy quark 
charged under U{1)d (see Section l2T2l) . The shaded regions in this figure are the same as in 
Figure [TH 

A particularly simple example is the one-flavor model, in which the A' is significantly 
lighter than the dark mesons and the physical dark Higgs boson associated with U{1)d 
breaking, hj:,. Under this assumption, the £"'"£"7 production is important, but off-shell pro- 
duction of dark-sector particles gives rise to high-multiplicity final states consisting of leptons 
and/or Standard Model hadrons that may yield more sensitive limits. 

For light A' , all dark mesons and ho can decay promptly to several y4"s. For example, 
the rj'jy decays through the chiral anomaly to 2A' . We show the 77^ decay length in its rest- 
frame in Figure [m for two different values of e, and in Figure [15] for the case that a^e^ is 
normalized to its DAMA/LIBRA expected value. Heavier dark mesons may decay to two or 
more A"s, or to an A' and a lighter dark meson. Baryon number is an accidental low-energy 
symmetry of the dark sector, so we will assume that the A/? is exactly stable, and all other 
dark baryons decay to it and one or more A"s. 

Therefore, all dark-sector final states eventually decay to a collection of y4"s and stable 
A/j-baryons. The A' usually decays promptly through kinetic mixing (cf. equations and 

my 

- -(^' ^ S X 10- c„, (l^) ' (1^) (^) , (36) 
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Figure 18: Top: Diagram showing 3-body decays of a heavy dark-sector excited spin-0 and 
spin-2 meson X£, into Standard Model states (through an off-shell A') and a lighter meson. 
Bottom: Decay length (cr) of Xd. Here an = o:, and e = 10"^ (left) or e = 10~^ (right). 
In the very light blue region (lower right triangle) the 2-body decay of Xd dominates. The 
remaining shaded regions in this figure are the same as in Figure [TTJ 



We expect to produce final states with at least one to two A' per dark hadron. These decays 
can easily generate high lepton multiplicities, even when relatively few dark hadrons are 
produced. 

In one-flavor models at high multiplicities, production of dark baryons is unavoidable and 
leads to invisible products which escape the detector. For low- multiplicity final states, when 
-\/i/m.^/ ~ 2, the efficiency for producing individual exclusive final states, such as those that 
do not contain any baryons, are reasonably high. 

Similar comments apply for multi-flavor models, in which the pions are significantly heav- 
ier than the A'. In this case, some pions may be stabilized by unbroken dark-sector flavor 
symmetries, or by an unbroken discrete subgroup of U{1)d if their charges are non- integer 
multiples of the dark Higgs charge. Kinetic mixing does not change the stability of these 
pions, and they will escape the detector. Pions that are not stabilized, however, can decay to 
2A' through the chiral anomaly, and these A'^s are observable through their prompt decays 
to Standard Model matter. 
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Figure 19: Decay length (cr) of a heavy dark-sector spin-0 and spin-2 excited meson X^, 
which has a 3-body decay to a hghter meson and to two Standard Model leptons or hadrons 
through an off-shell A' . Here = a, and the combination ao^^ has been normalized to 
its DAMA/LIBRA expected value, assuming that the dark matter is a 1 TeV dark meson 
neutral under U{1)ei, but consists of a light and a heavy quark charged under ?7(1)_d (see 
Section 12.21) . In the very light blue region (lower right triangle) the 2-body decay of X£, 
dominates. The remaining shaded regions in this figure are the same as in Figure [TTl 



3.2.4 Dark Sectors with Towers of Metastable Mesons and Baryons 

We now consider the one- flavor model with a heavy A . If the A' is heavy, then the tower of 
mesons lighter than 2m^/^ becomes metastable, i.e. they have no available decays within the 
dark sector. These decays are therefore suppressed by at least e^, and can often be displaced 
or even invisible. Likewise, all baryons lighter than m^^ + m^/^ are metastable for heavy A' . 

The discussion of dark pion decays in Section 13.2.21 applies to the r/^ as well. In partic- 
ular, decays must proceed through two A' gauge bosons, and if both are off-shell they are 
suppressed by e^. Therefore, they can easily be invisible. The decay of rj'j^ through one or 
two off-shell y4"s is included in Figures [T3] and [T3 

The decays of the heavier mesons are generically much easier to see. The vector meson 
generically mixes with the A', with a mixing angle that can be estimated from their masses, 

4 

m , 

~ , (37) 

This leads to a decay width that is suppressed by e^, and up to an unknown 0{1) coefficient 
is given by 




T(ujD ^ rn ~ — 1 _ ^ 1 + ^ . (33) 
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The total ujd width depends only mildly on m^^. In Figure [T6| we show cr for uoo for two 
different values of e, while in Figure [T7t we normalize a^e^ to its DAMA/LIBRA expected 
value. 

Exotic and excited mesons can have more spectacular decays. Spin-1 mesons are expected 
to decay into Standard Model pairs, just like the ujdi while spin-0 and spin-2 mesons can have 
3-body decays into Standard Model states (through an off-shell A) and a lighter meson, as 
shown in Figure [181 These have a typical lifetime, estimated by dimensional analysis, given 
by 



r(x, ^ A* + - + {uM) - i - (39) 

which can easily be displaced. This is shown in Figure [18] for two different values of e, and 
in Figure [19] for a case in which ao^'^ has been normalized to its DAMA/LIBRA expected 
value. 

Note that, unlike the decay of the r/^ into Standard Model states, where the 3-body 
region of phase space required a tuning of independent masses m^^ and m^/, the condition 
"^ri'^lujD < < 2m^/^ is satisfied naturally for a tower of exotic mesons Xj:,. Similar decays 
are expected for excited baryons. 

To summarize: although significant missing energy from both baryons and 's is gener- 
ically expected, typical events in such a sector also have a combination of prompt and dis- 
placed decays. 



3.3 Dark Sectors with Uncolored Fermions 

We have focused on Higgsed and confined sectors with no light fermions. However, the 
presence of a light fermion changes the phenomenology of either a Higgsed or a confined 
dark sector dramatically. Some or all of the bosons that are produced can decay to two dark- 
sector fermions, instead of decaying to Standard Model fermions through kinetic mixing. The 
dark-sector fermions are singlets under any confining gauge group (in the confined case, the 
dark-sector gauge group may have a subgroup which does not confine). The lightest dark- 
sector fermion does not have any decay modes induced by kinetic mixing, but can decay 
through new mixing operators of the form 

^{Li,h^^){LhsM), (40) 

where h^^ is a collection of dark-sector Higgses with the same quantum numbers as a dark- 
sector fermion L^, L is a Standard Model fermion, and hsM is a Standard Model Higgs. 
The decay lifetime of the dark fermions is quite model-dependent, but they will typically be 
long-lived. In this case, any dark-sector bosons that can decay to a pair of light dark-sector 
fermions may well be invisible. The only way to find such invisible sectors is by observing 
their recoil against other particles in the final state (in the context of B- factories, a photon). 

In dark sectors with more than one light fermion, this situation is ameliorated. The 
second-lightest fermion /2 may have decays, through either Higgs quartic or gauge boson 
kinetic mixing, 

/2^/iA'^/i£+r, (41) 
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where the A' may be on- or off-shell. This decay is especially important when the decay 
/2 3/1 is kinematically not allowed. These fermions generically have longer decay lifetimes 
than their bosonic counterparts by the ratio of three-body to two-body phase space, ~ 
l/lGvr^. Such decays are much harder to find than the bosonic analogues because they contain 
missing energy, and the invariant mass need not reconstruct a resonance. Instead, the 
dilepton invariant mass distribution may have an endpoint. Moreover, if a sector contains 
both metastable fermions and metastable bosons, the bosonic decay products may reconstruct 
resonances. 

4 Searches at low-energy e^e~ colliders 

Experiments at high-intensity e^e~ colliders at low energies, such as the B-factories BaBar 
and BELLE, can open a direct window into the physics of a low-mass dark sector by observing 
production of dark-sector particles and their decays back into the Standard Model. Though 
striking, these final states are quite different from those of interest for precision measurements 
of the Standard Model, and could have gone unnoticed in existing B-factory datasets. A 
common feature of the dark-sector production processes is a high multiplicity of leptons, 
pairs of which reconstruct resonances. However, the precise multiplicity and energy of these 
leptons, and the presence of displaced vertices or missing energy, depend on the detailed 
structure of the dark sector. 

This variety makes inclusive searches desirable, but the primary concern is maintaining 
acceptance for a large fraction of signal events while removing Standard Model backgrounds. 
To this end, we identify several search strategies, with exclusive search requirements for events 
with large physics backgrounds and inclusive treatments for more striking signal events. Our 
discussion is mostly qualitative; a quantitative study is beyond the scope of this paper. Our 
focus is on BaBar and BELLE, which represent the current high-intensity frontier, but similar 
analyses with data from KLOE, CLEO-c, the upcoming BES-III [53], and elsewhere could 
be more sensitive, depending on the nature of the dark sector. 

We identify six regions of interest, where different search strategies are required: 

• 4£ (exclusive), reconstructing Ecm (also 4£ -|- 7) 

• 4£ (exclusive), with displaced dilepton vertices (also 4£ -|- 7) 

• > 5£ + tracks (inclusive) , reconstructing Ecm (also + 7) 

• > 5£ + tracks (inclusive) , with displaced dilepton vertices (also + 7) 

• Very high track multiplicity, with many tracks consistent with leptons 

• 7 + nothing 

Before discussing each of these searches in turn, we list three prominent physics backgrounds 
that contribute to the channels with four or more leptons: 

1. QED processes, e+e^ Ai, which can reconstruct the full beam energy. 
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2. Two photons, emitted from the initial-state electron and positron, can combine to 
produce e~^e~, or tt+tt^ pairs. The electron and positron are usually, but not 
always, lost along the beam pipe. A 4-lepton event is obtained if the initial state pair 
is observed, or if an additional 7* produces two more leptons. The latter kind of event 
can be rejected by requiring reconstruction of the full beam energy. 

3. Sequential leptonic B decays produce events with up to four leptons, plus additional 
tracks; some energy is carried by neutrinos so these events do not reconstruct the full 
beam energy. 

4i exclusive, reconstructing E^ra (with or without 7) : Four-lepton final states {^Ji^^Ji~ 
/i'^/i", yU+/i^e+e^, and e+e'e+e^) that reconstruct the full beam energy can be produced 
by the simplest decay chains in Higgsed sectors (Figure [T3l(a)). or by confined sectors with 
m^^ < These processes also produce r and hadronic final states, but these are harder 

to search for, whereas muons have the lowest backgrounds. Requiring the reconstruction of 
the full beam energy Ecm removes the physics backgrounds 2) and 3) mentioned above, but 
also makes these searches insensitive to production modes that include stable or long-lived 
final states, such as baryons or metastable Higgs bosons. Pairs of leptons should reconstruct 
very narrow (resolution-limited) resonances, which can be used as a very useful discriminator 
between signal and background. The two lepton pairs need not reconstruct the same mass, 
but it is likely that they will reconstruct indistinguishable masses in an 0{1) fraction of 
events. Similarly, 4£ + 7 final states can be produced by simple Higgsed-sector decays, or by 
confined sectors with m^^ < rnA'/^ and itia' < Ecm- In this case, in addition to pairwise 
mass reconstruction as in the pure 4i case, the four leptons should now also reconstruct a 
(potentially broader) A' resonance. 

4£ exclusive with displaced vertices (with or without 7): The requirement of dis- 
placed vertices allows significant reduction of Standard Model backgrounds 1) and 2), without 
demanding that the observed products reconstruct the full beam energy -Ecm- In an exclu- 
sive four-lepton final state, background 3) is also negligible. Dropping the Ecm requirement 
is quite important — in scenarios with displaced vertices, it is also generic to have miss- 
ing energy coming from long-lived particles decaying outside the detector, or exactly stable 
baryons. For example, a confined sector with one light flavor can have displaced decays, 
and long-lived invisible 77^ decays. Similarly, the beam-energy requirement can be relaxed 
in a 4£ + 7 search and replaced by a displaced vertex cut. In this case, a four-lepton recon- 
structed mass peak is no longer expected as some of the A' mass can go into invisible decays. 
Pairs of same-flavor leptons should still reconstruct narrow resonances. 

> 5i exclusive (with Ecm or displaced vertices, with or without 7): More complex 
Higgsed-sector decay chains such as Figure [T3l(a). where the Higgs decays within the detector, 
or confined-sector decays with larger ratio Ecm/ can generate moderate track multiplicities. 
We stress that, though the leptonic final states are most distinctive, especially at higher 
masses the dark resonances will decay to r and hadronic states as well. Therefore, it is 
essential not to veto on non-leptonic tracks, and to tag leptons as loosely as possible to 
suppress backgrounds. It may also be interesting to specifically search for events with r+r" 
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or TT+TT^ pairs. Therefore, a more inclusive approach is warranted, i.e. requiring 5 or 6 leptons 
to reduce background, plus additional tracks. As in the M searches, two complementary cuts 
are useful for discriminating from background: either full reconstruction of Ecm or displaced 
vertices. 

High Track Multiplicities A confined dark sector with m.^^ ~ ^n'QCD Produces a high 
multiplicity ~ 10 — 20 dark hadrons, of which a significant fraction decay visibly within the 
detector. The resulting dark-sector parton shower can lead to somewhat coUimated jets of 
decay products. A large fraction (between 4/5 and 1/3 depending on meson mass) of dark 
mesons go into e or /i pairs, and this can be used to assist in distinguishing these final states 
from backgrounds. There is no source of analogous events from a Higgsed dark sector without 
a confined group. 

7 + nothing: A further interesting region of parameter space is when all the produced 
dark-sector states decay outside the detector, in which case all dark-sector products are 
invisible. This occurs, for example, if all states in the dark sector can decay to a long-lived 
pion or dark Higgs. In this case, the '7 + nothing' mode is the only effective search at 
low-energy e^e~ colliders (see also [17J). These events can be detected by observing a peak 
in the photon energy spectrum at a/s ~ m^'. 

The searches mentioned above are complementary to the dilepton resonance searches 
(e+e~ 'J A' 'yi~^£~) that have recently been discussed in [TTl [35l HI]. However, the 
A' decay is significant only when the A' is the lightest state in the dark sector. Even 

if the A' is the lightest state of the dark sector, any other kinematically accessible states are 
pair-produced with a comparable rate through an off-shell A' , which produces a signal in the 
multi-lepton channel, which has much lower QED backgrounds. The multi-lepton searches 
considered here overlap with the 4- and 6-lepton searches discussed in [16] in the context of 
Abelian dark sectors, but we emphasize that much higher multiplicities are easily obtained 
in more general non- Abelian dark sectors. 

5 Conclusions 

Evidence for a low-mass dark sector may be hidden in unexamined e~^e~ collider data. If 
the new sector contains a f/(l)z) factor, then moderately large kinetic mixing e ~ 10~^ is 
not only generated by generic mechanisms, but suggested by a variety of experiments and 
observations that look for dark matter. States of a new sector with couplings in this range 
can be efficiently produced in e~^e~ colliders, and may decay spectacularly, so that potential 
discoveries are within the reach of searches in existing data from low-energy e'^e~ collider 
experiments such as BaBar, BELLE, CLEO-c, and KLOE. 

The evidence for a "dark sector" from terrestrial and satellite searches for dark matter 
is indirect, but strikingly consistent. The large local electron/positron excesses reported 
by HEAT, PAMELA, PPB-BETS, and ATIC suggest interactions of 0(TeV)-mass dark 
matter with a new boson with mass 0{1 GeV). Simultaneously, two experiments suggest 
multiple dark matter states with small mass splittings: the INTEGRAL 511 keV-line and 
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DAMA/LIBRA modulation signal are consistent with 1 — 10 MeV and ~ 100 keV splittings 

among dark matter states, respectively. New forces, Higgsed or confined at 0{GcV), can give 
rise to these splittings. Any subset of these results suggests a departure from the standard 
minimal WIMP paradigm, but dark matter data will provide at best an indirect probe of 
such structure. 

If the DAMA/LIBRA modulation signal comes from inelastic scattering mediated by 
kinetic mixing, then it is the first quantitative probe of kinetic mixing between the Standard 
Model and the dark sector, and can be used to estimate the production rate of new dark- 
sector particles. If the mass of the new U{1)d gauge boson is uia' — 1 GeV, cross-sections 
~ 0.1 fb are expected in a Higgsed dark sector with an — oi. This estimate leads to 50-100 
events expected in the stored datasets from BaBar and BELLE. The e+e~ production cross- 
section normalized to DAMA/LIBRA scales as m\,, so much larger cross-sections ((9(pb)) 
are possible for a heavier A', while higher-luminosity experiments are required to observe 
dark sectors with lower niA' ~ lOO's of MeV. Cross-sections up to 100 fb are possible even 
for a 1 GeV A', if the DAMA/LIBRA scattering is suppressed by a form factor, as expected 
in confined models. 

Though spectacular, the signatures of dark-sector production at e+e~ colliders would still 
have been missed by existing searches. The precise final states depend on the spectroscopy 
of the dark sector — general regions of parameter space allow a combination of prompt and 
displaced decays, and states that escape the detector invisibly before decaying. All decay 
modes are expected to produce lepton-rich final states. Moderate multiphcities (4-10) of 
leptons and very high multiplicities are both possible, and require different search strategies. 
Some of these lepton pairs are expected to reconstruct very narrow resonances, which in 
principle can be combined to reconstruct the spectrum of the dark sector. Dark sectors with 
very long-lived final states may only be visible in a '7 -|- nothing' search. Further work to 
optimize the strategies proposed here and study their efficiencies is an important task. 

The small production cross-sections for uia' ^ 100 MeV pose a particular challenge. 
Novel experiments are called for to explore this region. Two features of the low-mass region 
are noteworthy. Lower-energy 0(100 MeV) colliders may reach higher luminosities than can 
be attained at GeV-energy coUiders. Moreover, not far below the cross-section reach of B- 
factories, at e ~ 10~^, the expected decay lengths of spin-1 states in the dark sector become 
0{1 m), so that very different techniques, such as beam-dump experiments, may be the most 
fruitful way to isolate dark-sector production events. 

Dark matter experiments suggest new low-energy gauge interactions beyond the Standard 
Model. If a dark sector exists, it dramatically reshapes our understanding of the structure 
of nature. Searches in existing collider data offer a simple but powerful probe of this new 
dynamics and should be pursued. 

Acknowledgments 

We are grateful to Daniele Alves, Nima Arkani-Hamed, Siavosh Behbahani, Howard Haber, 
Jay Wacker, and Neal Weiner for many fruitful discussions. We especially thank Mathew 
Graham and Aaron Roodman for numerous discussions regarding BaBar physics, and Michael 
Peskin for helpful feedback on our analysis of production and decay phenomenology. RE and 
PCS are supported by the US DOE under contract number DE-AC02-76SF00515. 



35 



A Constraints on Kinetic-mixing between hypercharge 

and a new U{1)d 

In this appendix, we review the most important existing constraints on the kinetic mixing 
between Standard Model hypercharge and U{1)d. 

For low-mass (MeV-GeV) dark photons A', there are a variety of constraints on the 
kinetic mixing that are summarized in [35] (for a discussion of searches for even lower mass 
dark photons see [5ll|55]). Constraints from Big-Bang Nucleosynthesis disfavor tjia' < 0(few 
MeV) [56]. For ttia' ^ 10 MeV, the strongest constraints come from the A' contribution to 
the anomalous magnetic moment of the electron, , while for m^/ above 10 MeV, but less 
than a few GeV, the tightest constraints come from the ^'-contribution to the anomalous 
magnetic moment of the muon, . For m^/ ~ (9(100 MeV - 300 MeV), constraints from 
certain rare decays of mesons, such as tt^A', are in some cases competitive with the 

constraint on , but depend sensitively on the decay modes of the A' — see [22] for a 
more detailed discussion. Above a few GeV, kinetic mixing is only very weakly constrained. 
However, since kinetic mixing also induces a coupling between the Z and dark-sector matter, 
a potential probe could come from searches for rare Z-decays. 

Since we are interested in the mass range m^/ ~ 10 MeV - 20 GeV, we focus on the 
constraint on and the potential probe coming from rare Z-decays. 

Anomalous Lepton Magnetic Moments: Let us first discuss the contribution from the 
A' to the anomalous magnetic moment of a lepton, i. This was calculated in [35], and is 
given by 

a? = / dz ^, ^ \, . (42) 

For a^', this can be converted to a constraint of < 10~^(m^'/(10 MeV))^ [35], which is 
important for m^/ < 10 MeV. A constraint on is more ambiguous, since the theoretical 
prediction for the muon anomalous magnetic moment within the Standard Model, a^, involves 
a hadronic contribution which must be estimated from experiments, which do not all agree. 
Using data from e^e~ annihilation to hadrons, the theoretical value of is smaller than the 
measured value by (302±88) x 10^^^, a ?>Aa deficit [23 EH]- However, estimates involving data 
from r's or from preliminary BaBar results on the precision measurement of the e^e~ 
n^Ti~{^) cross-section with the ISR method |i58] do not point to a significant discrepancy. 
Since the A' contribution to is positive, a conservative constraints on the e— m^' parameter 
space is obtained by taking the estimate from e^e~ annihilation to hadrons. In particular, 
for the constraints shown in Section [21 we follow [2S] and require < (302 + 5cr) x 10^^^ = 

7.4 X 10^^. For m^/ ^ m^, equation reduces to ~ a^^^' ^"^^ constraint on 
thus becomes very weak for m^' larger than a few GeV. 



Exotic Z-Decays: For ruA' larger than m^, searches for exotic Z-decays among the 2 x 10'' 
Z-bosons recorded by the experiments at LEPI could potentially probe kinetic mixings more 
sensitively than a^. Let us discuss this in more detail. The couphng induced between the Z 
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and dark-sector fermions is given by 

Nf 

Cd-J2^9d tan OwZ^.X-f'^X, (43) 
1=1 

where Nf is the number of dark-sector particles couphng to A' with gauge couphng qd, and 
sin^VK is the Standard Model weak-mixing angle. In this expression, the A' coupling to 
the electromagnetic current has been normalized as £ D eg^A'^J^j^, see equation ([1]). The 
Z-decay width to dark-sector fermions is thus 

r{Z XX) ~ l-N^e'anmz t^n'ey, V Jl-^(l + (44) 
3 ^ y m| V m| / 

where Nc is the number of colors of the dark-sector gauge group. 

Since the total width of the Z has been measured to be 2.4952 ± 0.0023 GeV [39], the 
branching ratio oi Z XX is given by 

BR{Z ^ XX) ^ 2.8x10-^ N,— (-^)\ (45) 

a VIO 

where we assumed mx rriz- Depending on the decay modes of X, which we discussed 
in detail in Section [3l a variety of signatures are possible. For example, if the X's decay 
invisibly or outside the detector, the constraint is not better than the uncertainty on the 
width of the Z, which is about 0.1%. More exotic Z decays into multi-leptons are possible, 
but require dedicated searches. The best constraints on rare Z decay branching ratios are 
no better than the 10~^ — 10~^ level [39], and even dedicated searches may not reach this 
level for complicated final states. Equation fH5l) shows that for typical parameters that give 
0{fb) cross-sections at B-factories, the branching ratio Z —>■ XX is about three orders of 
magnitude smaller than what could be probed at LEPI even with a dedicated search. For 
the multi-lepton final states we consider in this paper, the B-factories will thus be a far more 
sensitive probe of the e — m^/ parameter space than LEPI, at least for rriy^/ less than a few 
lO's of GeV. 



B DAMA/LIBRA Normalization of Couplings 

Inelastic dark matter offers a simple explanation that reconciles the annual modulation signal 
reported by DAMA/LIBRA with the null results of other experiments, and is reviewed in [T^ 
[T3j . Below, we summarize our procedure for normalizing aoe"^ assuming that DAMA/LIBRA 
is explained by an iDM mechanism. We will normalize a^e^ for two iDM scenarios: one in 
which the inelastic channel proceeds through charged-current scattering as in [12], and one 
where the inelastic channel is through a hyperfine transition as in [9J . 

For charged-current scenarios, the inelastic cross-section to scatter off of a nucleus of mass 
m^v and charge Z recoiling with energy Ept is 

da _ SnZ^aane^m, 1 ^^^^ 



dEn {2mNER + m\ 



A' 



37 



where a is the fine structure constant, v is the dark matter nucleus relative velocity, and m^/ 
is the A' mass. F{E^ is the Helm nuclear form factor [59], and it accounts for the loss of 
coherence scattering off of the entire nucleus at large recoil. We use. 



|F(i5H)P=(^^f^ye-^^l^^ (47) 

where s = 1 fm, Tq = — 55^, r = 1.2 A^^^ fm, and \q\ = y/2mj\fE£i. 
For hyperfine transition scattering |9], the inelastic cross-section is 



da ^ 47rZ^aCi„Q;De^ mj^ER 
dEn ^ v^Al {2niNER + ni\ 



\F{Er)\\ (4J 



where Cmj ~ GeV ^ is the N electric dipole of the dark matter. To fix the relation 
between and the hyperfine splitting 5 in this case, we introduce a scale A^j such that 

^ - ^ (49) 



As discussed in Section I2l2l K^)! K^) is naturally a small number. We will therefore normalize 
a^e^ as a function of m^/ for different choices of the small dimensionless coefficient, 

C=|^l. (50) 

Following [60], the differential scattering rate per unit detector mass is, 

- ^ d\ -fiv;v.)4^, (51) 



dEji mNniDM Jv^uer) dE^ 

where p = 0.3 GeV/cm^ is the local density of dark matter, and f{v) is the dark matter 
velocity and velocity distribution function in the lab frame. Introducing a 0(fesc — v) to 
naively cut off the velocity profile above the galactic escape velocity Vesc ~ 500 — 600 km/s, 
we use a simple velocity distribution function. 



f{v;ve)(x[e -e ^'o jQ{vesc - \v - Ve\), (52) 

where vq ~ 220 km/s, and Ve is the earth's speed in the galactic frame. We use the numerical 
values reported in [60] for Ve. For the annual modulation rate, we use 

(53) 



/ dR ' 




^ dR 




dR 


[dEn, 


modulated 


KdEjt 


June 2 


dER 



Dec 2 , 

To normalize the combination a^e^, we fit to the reported DAMA/LIBRA signal, using the 
combined 0.82 ton-yr exposure |10i|. We do this by minimizing a function using the 12 
half-keVee bins and reported error bars between 2 — 8 keVee with 10 = 12 — 2 independent 
degrees of freedom. The lower and upper ranges of a^e^ are taken by scanning over the 
region with ^ 16, and typically differ by a factor of ~ 10. Our estimate is obtained by 
taking the geometric mean of the lower and upper values of a/^e^. 
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